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We examine the dependence of the physical quantities of the rotatory molecular motor, such as 
the rotation velocity and the proton translocation rate, on the chemical reaction rate using the 
model based only on diffusion process. A peculiar behavior of proton translocation is found and the 
energy transduction efficiency of the motor protein is enhanced by this behavior. We give a natural 
explanation that this behavior is universal when certain inequalities between chemical reaction rates 
hold. That may give a clue to examine whether the motion of the molecular motor is dominated by 
diffusion process or not. 

PACS numbers: 05.40.-a, 05.60.-k, 87.15.-v 

I. INTRODUCTION 

Many kinds of molecular motors are known to play essential roles of life. The word 'molecular motor' is, in a 
wide sense, used to describe protein molecules that transduce chemical energy from some source into mechanical work 
(and in some case, vice versa). Molecular motors are classified into two kinds: linear and rotatory motors. As the 
former, myosins, kinesins, and dyneins, which slide along the specific filaments in the specific direction, respectively, 
are known to work for muscle construction, transport of materials in cell, etc. by hydrolysing ATP 0,0- As the latter, 
two motors are known to date. The one is bacteria flagellar motor, which works as a propeller by gaining chemical 
energy from the proton gradient across the membrane. The other is ATP synthase which is found in membranes of 
living things. It works wholly as a transducer between the energy of proton gradient across the membrane and that of 
ATP synthesis from ADP and inorganic phosphate. It can both synthesize ATP and pump proton against the proton 
gradient by hydrolysing ATP (see Fig. 16-28 in Ref.|3|). It is composed of two parts, the F a - and the Fi-part, each is 
known to work as a rotatory motor 

The F -part is embedded in the membrane. It contains a proton channel 
and transduces the energy of transmembrane proton gradient into rotatory torque when working as synthesizer. The 
Fx- part plays a role of ATP synthesis/hydrolysis, transducing rotatory torque into the energy required for ATP 
synthesis and vice versa. 

It is believed generally that the energy transduction efficiency of molecular motor, that is, the ratio of the energy 
output to input from the energy source, is extremely high. For example, it is reported in Ref.|l2| that the efficiency 
is almost 100 percent for the Fi-part of ATP synthase. Thus, understanding the mechanism of molecular motors, 
especially that of energy transduction is important not only because of their many roles of life, but also because of 
the possibility that, if this mechanism is novel, it may be applied to a new technology as a machine with high energy 
transduction efficiency, even if it can be used only in a constrained scale. In the scale of motor molecule(^lO-lOOnm), 
water molecules collides motor molecule so many times. These collisions make the diffusive motion of the motor 
molecule and thus it is natural to infer that this collision effect contributes to the motion of motor, but there is no 
crucial evidence that it is essential. 

Here we take up a simple model inspired by the F^-part of ATP synthase|0, |2(J and examine the dependence of 
the rotation velocity and the proton translocation rate of this model on the chemical reaction rates. The efficiency of 
energy transduction is calculated from these quantities. As far as we know, there is no computation of this kind. 

II. MODEL 

The -Fo-part is schematically drawn in Fig's^and^ It contains three types of subunits, a, b, and c. As shown in 
FigH the a-subunit is fixed in the membrane. The b-subunit connects the F a - and Fi-part but is not shown in Figs^ 
and[2Jsee Rcf. 3j). The c-subunits are arranged in a ring and this ring rotates accompanying the proton translocation. 
There is a proton binding site in the almost middle of each c- subunit, Asp61, a carboxyrate. A proton passes through 
the interface between the a- and c- subunits (but as described later, not directly), called proton channel. In the proton 
channel, there are two paths to these binding sites, the left one is from the basic side, the inner side of the membrane, 
and the right is from the acidic side, the outer(see Fig|2J)[3- Proton concentration in the acidic side are kept higher 
than that of the basic side by respiratory chains, which pumps protons in the basic side out of the membrane. Roughly 
speaking, a proton flows into and binds to the right site from acidic side, goes through the membrane with c-ring 
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FIG. 1: The proton channel viewed from the basic side, the inner side of the membrane. There are the c-subunits arranged in 
a ring and the ring rotates counterclockwise during synthesis. The position variable x is defined as a rotation angle of the ring. 



FIG. 2: The proton channel viewed from the a-subunit. There are two proton binding sites in the channel and only these two 
can be either protonated or unprotonated. The left path of protons is from the basic side and the right path is from the acidic 
sideQ. These proton binding sites move rightward during synthesis despite the leftward load from the Fi-part. 



rotation, and dissociates from the left site and flows out to the basic side. That is, a proton passes through the proton 
channel not directly but via the membrane accompanying the rotation of the c-rin g. 

We investigate the above situation using the '(simply) biased diffusion model' jlij. l2*cj. It was first presented to 
explain the rotation of bacteria flagellar motor which has similar structure to the i^-part of the ATP synthase but is 
one order larger in linear dimension. The mechanism is as follows: There are N proton binding sites in the rotor (i.e. 
the ring is composed of N c-subunits and iV=9-12 in Escherichia coli) at same intervals and two of them are in the 
proton channel. Only two proton binding sites are allowed to be in the channel and each of them can be protonated 
or unprotonated. The others can be only protonated due to the hydrophobicity of the membrane. Thus the states 
of the motor are determined by the protonation of the two proton binding sites in the channel: E is the empty state 
where both sites are unprotonated, R is the right state where the right site is protonated and the left is unprotonated, 
L is the left state where the left site is protonated and the right is unprotonated, and F is the full state where both 
sites are protonated. We assume that a proton does not hop between two proton binding sites. When one of the 
proton binding sites in the channel comes to the boundary, it can move into the membrane only if it is protonated 
and otherwise cannot due to hydrophobicity of the membrane. If there is difference from detailed balance between the 
rates of one site at which a proton binds to or dissociates from it and those of the other due to the transmembrane 
proton gradient, it become possible to take out a net motion. 
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Dynamics of each state is described by the Langevin equation, 



dx 
~dt 



-7r 



d (j)i{x) 
dx 



Ri(t), 



(1) 



where i refers to the state, 7$ is the friction constant, x is the position of the motor (see Fig.l) which is explicitly 
defined later. The potential 4>i(x) represents the load torque r from the F\ part, so (f>i(x) = tx \l2\. Ri(t) is the 
random force of Gaussian noise, which satisfies 



(Ri(t)) = 0, 



(2) 



(R i {t)R j (t')) = 2D i 5 ij S(t-t'), 



(3) 



where (• • • ) is the time average and Di is the diffusion constant satisfying the Einstein's relation, Di = jik^T 
(T: temperature, k^: the Boltzmann constant). The constants Di and 7^ depend on the states since the different 
conformations have different surface areas, densities, etc. But for simplicity, it is assumed below that all the states 
have effectively the same values, 7, = 7, and Di — D, since those differences are small. In addition to the above 
single-state motions, transitions between different states occur. Thus it is convenient to adopt the formulation of the 
Fokker-Planck equation equivalent to Eq. , 



di 



p(x,t) 



d_ 

dx 



Tl(x,t)+K(x)-p(x,t), 



(4) 



U(x, t) = 7f (x)p(x, t) - D—p(x, t), 



(5) 



where f is the external force f = diag[r, r, • • ■ , r]. K is the transition matrix which describes changes between states. 
The system discussed here has 4-component probability p(x,t) and flow Tl(x,t) 



p(x,t) 



p E (x,t) 

Ph{x,t) 
p F (x,t) 



n(x,t) 



tt E (x,t) 
U R {x,t) 

n L {x,t) 

U F (x,t) 



(6) 



where Pi(x, t) and Tii{x, i)(i=E,R,L,F) describe the probability and its flow that the motor in state i is at position x 
and at time t. The coordinate x is placed as follows: The origin O is set at the center of the c-ring and the position 
variable x is defined as a rotation angle (Fignj. The proton channel extends from x — to i5 (5 — 2tt/N and N = 12 
fixed hereafter) . The position x — is defined that the left proton binding site is at the left boundary of the channel, 
and x — 8 is that the right site is at the right boundary. When one site moves into the membrane, another site appears 
into the channel from the other side of the membrane. In this sense, the system has periodicity. These constraints of 
periodicity and hydrophobicity mentioned above are expressed by giving Eq.Q the boundary conditions, 



n E ((M) = n E (5,t) = 0, 



(7) 



and, 



The transition rate matrix is written as 
K = 



n R (o,t) = n L ((5,i) = o, 
n L (o,t) = n R ((5,i) ) 
n F (o,i) = n F (<M). 
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(8) 
(9) 
(10) 

(11) 
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There are a few assumptions in the above expression of the transition rate matrix : 1) The chemical reaction, that 
is, binding and/or dissociation of proton to/from the sites, is sufficiently fast compared with the motion of the motor 
protein. This is justified because the mass and size of a proton is much smaller than those of the motor so the diffusion 
coefficient of a proton is much larger than that of the motor protein. 2) There is no correlation between the reaction 
of the left site and that of the right site. 3) The reaction rates are independent of the position of the motor x. This 
assumption may have to be modified to take account of other influences, for instance, interaction between residues|l9j|. 
But we neglect the influences like these for simplicity. In this case, for example, it is thought that &er, the rate at 
which state E switches to state R and feu, the one a t which L switches to F are both written by kj^, the one at which 
a proton binds to the empty right site. So similarly we set, 

(12) 
(13) 
(14) 



&ER 


= &LF 


- k R 
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and, 



fcr.E = &fr = k 



(15) 



For the motor to work, the detailed balance should be violated due to the free energy acquired from a proton passage 
of the channel, AG. Thus the ratio of transition rate from L to R via E or F to that from R to L is written as, 



(L — » R) _ /CLEVER + &LF&FR 



(R — > L) /crenel 



FL 
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= exp[AG/fc B T], 



(16) 



where AG is a function of ApH (ApH = pH-Q — pHA, A, B denotes the acidic and basic side, respectively) and V, the 
membrane potential accompanying this difference of proton concentration and written as, 



Under these assumptions, 
the matrix Q, 



AG = V + k B Tln[10 ApH ]. (17) 
this model can be solved analytically. The transition rate matrix K is diagonalized by 
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(18) 



Q 



1 


1 


1 


1 


-k L 


-k L 

K in 


k L 

"■out 


k L 

""out 


uR 

ft in 


k R 
"out 


k R 

ft in 


k R 

"-out 


k L k R 
"In "'in 


jl,L L.R 
"'in ""out 


hit jLR 
"-out ""in 


jLL lR 
"-out "-out 



(19) 



and 



Q ^KQ 



-(kf n + k^ ut ) 



Then in a steady state, Eq.(@J is reduced to 
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The solution of Eq.JSlJ is, 



where d (i — 1, 2, • 



are integral constants and, 



Ci + C 2 e«* 
C 3 e v £ x + C 4 e'>i> 
C b e< x + C 6 e^ x 
C 7 e^ x + C a e^ x 



7T 

5 D' 



(22) 



(23) 



± _ 7 r± V( 7 r)^+4£>(fc^ n + fc^) 
~ 2# 



(24) 



and, 



Therefore, we obtain 



±_ 7 t± V( 7 r) 2 +4J(fc^ + fcH~) 
^ 2# 



± _ IT ± V(7t) 2 + 4ggS + fcjt + + *&t) 



(25) 



(26) 



p(x) = Q(Q _1 p(a;)) 



fc?ut^out(^i + Cae^) - ^(Cae"?" + C A e^ x ) - k^ ut (C 5 e< x + C 6 e^ x ) + C 7 e^ x + Cge'W 
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(27) 



and, 



n(x) = Q 
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(28) 



The integral constants Cj's are determined to satisfy the boundary conditions lf?f- l|10(l . Five of these six boundary 
conditions are independent but the rest one is not because from Eq.Q, the following relation is always satisfied in 
the steady state, 

d f s 

= — J dx^p l {x,t) 



{Ui(0,t)-Ui(6,t)). 



(29) 



In addition to these boundary conditions, the periodicity of the probability at the boundaries, 

Pr(*)=Pl(0), 



(30) 
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PF (S) =pf(0), 



(31) 



and the normalization condition 



J dx^p l {x) = l 



(32) 



are necessary in order to close the equations for Cj's. Using these conditions CVs are determined as, 

c = M- 1 s, 

where 
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(36) 



We find useful formulas for the probability and the flow in state i. For the probalities, the following relations hold 
in the steady state, 
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(37) 
(38) 
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where 



Pi = / dxpi(x). 
Jo 



(41) 
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The relations above are obtained by integrating Eq.0} over x and using the boundary conditions JJJ-^JJJ. For the 
flow of probavility, the following relation holds, 



dxllj (x) 



jTfi(x) + D\pi(0) -Pi(5)], 



(42) 

which are obtained by integrating Eq.JSJ with respect to x. 

Let us now derive formulas for the observables; the average rotation velocity (v), the proton translocation rate 
N(H + ), and the efficiency of energy transduction. The average rotation velocity (v) is given by the total sum of the 
integrated flow of probability , 



(v) 



(43) 



i=E,R,L.F 



^out)- 



(44) 



Using Eq.JSSJ, we obtain 

The proton translocation rate 7V(H + ) is given as follows: A proton goes out from the left binding site to the cytoplasma 
when the state L becomes E or F becomes R. And when E becomes L or R becomes F, a proton comes back to the 
left binding site. So the flow of the left binding site Jl is defined as 

j l = fco U tPL - kf n p E + k^ ut p F - kf n p R . 



In a similar way, the flow of the right binding site is written as 

JR = -koutPR + k hiPE - fcoutPF + k ^p L . 

For the system to be in a steady state, Jl and Jr. must be equivalent and we obtain 
N{R+) = J L = J R 

= ^[(fc£ ~ k in)PE ~ (kout + kf n )pR + (fc£ + fcou t )PL + - kf ut )p F }. 

Using Eqs. (|3^|l and i|39"|) . we find that this A^(H + ) is simply given by 11^(5) and evaluated as, 
A^(H+) = i[Eq.C3-Eq.G9] 

= n R (5)(= n L (o)) 

= fc£fcoutCi7T - k? n (C 3 D^eri s + C 4 D V +eS s ) 

+kt t (C 5 D7 lR e< s + C 6 Dr,+e^ 5 ) - (C 7 D^ R e^ s + C s D V +e^ 5 ). 



Finally, we define the efficiency of energy transduction as follows, 

-t(v) 



(45) 
(46) 



(47) 



This definition coincides with that given in |18| . 



AG • iV(H+) 



III. RESULTS 



(48) 
(49) 



We investigate the dependence of the quantities, (v), N(K + ) and e on the transition rate K using the above 
analytical solution. Parameters used for calculation are, according to Ref.0], D = 2 x 10 4 rad 2 sec _1 , k^T =4pN-nm, 
7 = 5 x 10 rad sec pN nm and r = — 40pN-nm. The transition rates are given as 
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(50) 
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FIG. 3: The rotation velocity (v). K of the abscissa is that of Eg. (l5UII . The solid line with +'s is the plot for pK a —4.5 and the 
broken line with x's is for pK a =5.5. There is almost no rotation in the region of small K. The faster the chemical reaction, 
the more this quantity increases monotonically and it saturates finally. 
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K 

FIG. 4: The proton translocation rate JV(H + ). The solid line with +'s is the plot for p7^ a =4.5 and the broken line with x's is 
for pK a ~5.5. There is almost no translocation observed for small K in both cases. For pK a = 4.5, there exists a region where 
the increase becomes very slow (K -~10-12). At large K, it diverges as 10 K ^ 2 . 



where pHx = 7.0 and pH^ = 8.4 are the proton concentrations of each side (A:acidic, B:basic, respectively), pK a 
is the acidity of proton binding site, <f> — 2.3k^T is the surface effect, and V = b.6k^T is the membrane potential 
(the proton binding site is assumed to be in the middle of the membrane). The overall factor 10 K means the proton 
absorption rate of the path. Since it depends on the proton diffusion coefficient and the surface area of the channel, 
we assume that this factor is common value for all transition rates. Variation of K is also regarded as varying pH of 
each side and pK a , preserving the differences of each other. 

We plot the dependence of (v), N(R + ) and e on K in Figs|3-[SJ For the rotation velocity, (v), no qualitative difference 
dependent on pK a is found. There is almost no rotation in the desirable direction when the chemical reaction rate 
is low (namely when K is small). Then the velocity monotonically increases as K increases and saturates when 
the reaction is very fast. For the proton translocation rate, A r (H + ), a qualitative difference dependent on pK a is 
found though almost no translocation is observed during small K in both cases. For pK a =5.5, N(H + ) increases 
monotonically. On the other hand, for pK a —A.5, while it increases monotonically also, but there exists a region where 
the increase becomes very s\ow(K ^10-12). In both pK a 's, N(H + ) diverges as 10 K / 2 when K becomes large enough. 
The efficiency, e, reflects the above results of (v) and 7V(H + ). In both pK a 's, it increases as K increases at first. Then 
it reaches a peak value (at K ~12 for pK a =5.5 and at K ~13 for pK a =4.5) and decreases monotonically afterward. 
The efficiency goes to zero as K diverges reflecting the fact that (v) converges to the finite value but A^(H + ) diverges 
for large K. Note that the peak value of the efficiency for pK a = 4.5 is enhanced and is about 25 percent larger than 
that for pK a = 5.5. This comes from the slow increase of iV(H + ) for pK a = 4.5 described above. 

Inferring from the forms of Eas. l|23|l - H26|l . the key to account for these results are thought to be the inequalities 
between the chemical reaction rates 's(i=E,R,L,F and j=in,out) and the reciprocal of the relaxation time of this 
convection-diffusion system, (jt) 2 /D. The inequalities between fcj's for pK a =5.5 are different to those for pK a =4.5. 
For pK a =5.5, the inequalities are fc^t > fc £ > fc Sit > while for pK a =A.5,k^ ut > fc^ ut > > k^ n . Under the 
parameters used here, which are plausible in living body, only the above two sets of inequalities are possible for the 
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FIG. 5: The transduction efficiency e. The solid line with +'s is the plot for pK a =4.5 and the broken line with x's is for 
pK a =5.5. It has a peak at K ~13 and e ~0.5 for pK a =4:.5, and K ~12 and e ~0.4 for pK a =5.5. For large K, it converges to 
reflecting the fact that (v) converges and 7V(H + ) diverges. 



model to work well. Note that the inequalities k^ > k^ and k^ ut > k^ ut hold since the proton concentration in the 
acidic side are higher than that in the basic side and the difference makes the transmembrane electrostatic potential, 
and also that k^ ut > k^ is required for the motor not to diffuse leftward by the load torque r from the Fi-part. 

While the above two sets of inequalities between fc'-'s do not depend on K, but the inequalities between (7T) 2 / D 
and kj's do depend on K. And fc* 's which are larger than (jt) 2 /D are relevant to the dynamics of the system. When 
K is small enough, (jt) 2 /D is the largest among them. Thus all the transition rates between the states are neglected, 
so the motor does not work well. As K increases, k^'s also increase, and when K ?»9.5 for pK a =4.5 and K «10.5 for 
pK a =5.5, the largest rate k^ ut becomes comparable to (jt) 2 /D. A proton in the left channel can dissociate before 
diffusion and state R is created. Thus IIr increases and is the largest among IV s in this region of K. 

Next, for pK a =4.5, k^ ut becomes comparable to (7T) 2 /D when K wl2. Near this region of K, a proton in the right 
channel in state R can dissociate before diffusion, so state R becomes state E. Thus fin decreases until YIe becomes 
larger than fl^. As a result, the increase of 7V(H + ) (— IIr(5)) is suppressed in this region, but (v) does not because 
the increase of fl^ cancels out the decrease of ILr (See Figs |3] and On the other hand, for pK a =5.5, ft^ continues 
to increase since kj^ is larger than k^ ut . Thus N(H + ) and (v) increase monotonically. 

At last, when K >16, all the transition rates become greater enough than (7r) 2 /D, so the matrix M is reduced to 
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where 
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FIG. 6: The integrated flow of state E. The solid line with x's is the plot for pK a =4:.5 and the dotted line with +'s is for 
pK a =5.5. It increases monotonically both for pK a —4.5 and 5.5. 
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FIG. 7: The integrated flow of state R. The solid line with x's is the plot for pK a =4.5 and the dotted line with +'s is for 
pK a =5. 5. For pK a =5.5, it increases monotonically. But for pK a =4:.5, there is a region where it decreases. 



and 

fLR = yJW n + % ut + k? a + kn t )/D. (54) 

(In this limit, rjf -> ±rfc.) From this, it is easily found that Ci_ 2 = 0(1O~ 2K ), C 3 _ 6 = O(10~ K ), and C 7 _ 8 = 0(1), 
so (v) converges when K diverges. The convergent value of (v)k— >oo is given by 

7 2 t 2 (5[-&c + a{e~< TS l D + b(e^ s / D - 1)}] 

~ bc{D(ef rS / D - 1) - -/tS} + a[D - De^ r5 / D + jTS{e^ s / D + b{e^ rS / D - 1)}] ' 

(55) 

where a = kf n /k^ ut , b = fe^/fc^ ut , c = fc^ t /fc£ ut . Comparing (gHJ) with (JUJ), it is also found that N(H+) = O{10 K/2 ). 
The divergent behavior of 7V(H + ) is consistent with Fig. 0] 

We show the integrated flows of probability for state E and R, IIe and fin, in Figs and Their behaviors 
are consistent with the above explanations. We also see that LTp and II l are much smaller than II e and IIr under 
the above two sets of inequalities (k^ ut > fc? > k^ ut > k\ n and k^ ut > fc^ ut > A;?; > fc^ ). (See FigOO) Thus the 
qualitative behavior of (v) and N(H + ) is determined by only the states E and R. This is also consistent with the 
above explanation which neglects the state L and F. 

IV. SUMMARY AND DISCUSSION 

We have investigated the rotatory molecular motor using the simply biased diffusion model. The model depends 
on four chemical reaction rates, k™, k^ ut , k 1 ^, &^ ut , the diffusion and the friction constants of the motor, D, 7, and 
the load torque, r. We have solved the model analytically, and examined the relation between these chemical reaction 
rates and the physical quantities, such as the rotational velocity, the proton translocation rate, and the efficiency of 



11 




FIG. 8: The integrated flow of state L. The solid line with x's is the plot for pK a =4.5 and the dotted line with +'s is for 
pK a =5.5. It has a peak for K ~11 for pA" a =5.5 and K «10 for pK a =4:.5. 



the energy transduction. It is found that below such value K that k^ ut , the largest among fcj's, becomes comparable 
with (7T) 2 /D, both (v) and 7V(H + ) are very low so this model does not work well and result in a bad efficiency and 
that there exists an optimal value for the chemical reaction to maximize the efficiency of the energy transduction. We 
also found that when the inequalities k^ ut > (jt) 2 /D > k^ ut > k^ > kf n hold, the increase of the proton translocation 
can become very slow and the efficiency of the motor is enhanced. This efficiency enhanced mechanism is naturally 
explained from this set of inequalities. 

This model is based on the diffusion process like the thermal ratchet model for linear molecular motors j^J, so that 
the efficiency enhance mechanism above may give a clue to examine whether diffusion process dominates the motion 
of molecular motor essentially or not. The inequality is controlled by the dissociation constant of the proton binding 
sites, pK a , the proton concentration of both side of the membrane, the proton diffusion constant, etc. 

The efficiency of this model is still not high comparing with that reported in Ref. |l2j | . There are some possibilities 
which can explain this: 1) The definition of efficiency may be still vague and controversial. In Ref.yjh the authors 
adapted a different definition of the efficiency from ours. Our definition coincides with that in Ref.[lj|, where the 
energetics of the Langevin equations and the Fokker-Planck equations are carefully analyzed. 2) There may be a 
necessity to reconsider the several relations between physical quantities which we assume to be appropriate near 
equilibrium, such as the Einstein's relation, since molecular motor is a non-equilibrium system far from equilibrium. 
3) There may be another mechanism which enhance the efficiency. For instance, an electrostatic interaction between 
residues is discussed in Ref.|l9j. Furthermore, experimentally conformational change of the motor protein during 
motion is reported 22]. It is important to construct a model which gets along with this experimental result. These 
may be the key to account for the high efficiency of energy transduction of molecular motors. 
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